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ABSTRACT 


Presently,  the  United  States  Navy  is  searching  for  an 
improved  method  to  predict  the  damage  to  a  ship  hull  or 
underwater  structure  that  results  from  an  underwater 
explosion.  One  method  of  predicting  this  damage  is  through 
the  use  of  nonlinear  finite  and  boundary  element  analysis. 
Underwater  Shock  Analysis  (USA)  code  combined  with 
VEC/DYNA3D  code  is  used  for  the  analysis  of  the  effect  of 
explosive  shock  on  numerical  models.  Initial  geometric 
imperfections  are  introduced  in  the  numerical  model  using 
modal  imperfections.  The  resulting  numerical  model  is  then 
subjected  to  a  simulated  underwater  shock  using  the  combined 
USA/DYNA3D  code.  A  sensitivity  analysis  is  performed  to 
look  into  the  details  on  the  damage  resulting  from  these 
simulations . 


TABLE  OF  CONTENTS 


I .  INTRODUCTION  .  1 

II.  INITIAL  IMPERFECTIONS  .  2 

III.  TWO  DIMENSIONAL  INFINITE  CYLINDER  MODELS  .  5 

A.  Perfect  Cylinder  Model  .  7 

B.  Imperfect  Cylinder  Models  .  10 

IV.  THREE  DIMENSIONAL  RING  STIFFENED  INFINITE  CYLINDER 

MODELS . 18 

A.  Perfect  Cylinder  Models  .  18 

B.  Imperfect  Cylinder  Models  .  21 

C.  Parametric  Study  on  the  Effect  of  Cylinder 

Geometry . 2  9 

V.  THREE  DIMENSIONAL  RING  STIFFENED  FINITE  LENGTH 

CYLINDER . 51 

A.  Perfect  Cylinder  Models  .  53 

B.  Imperfect  Cylinder  Models  .  59 

C.  Exponential  Decay  Shock  Wave  .  72 

VI.  SUMMARY  AND  CONCLUSIONS . 81 

APPENDIX  A:  FORTRAN  PROGRAMS  FOR  MODIFYING  INGRIDO  FILE 

FOR  MODAL  IMPERFECTIONS  .  83 

APPENDIX  B:  INGRID  INPUT  FILE  FOR  TWO  DIMENSIONAL  INFINITE 
CYLINDER  MODEL  .  89 

APPENDIX  C:  MATERIAL  MODEL  DETAILS  .  90 

APPENDIX  D:  INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  RING 

STIFFENED  INFINITE  CYLINDER  MODEL  .  92 

APPENDIX  E:  INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  RING 

STIFFENED  FINITE  CYLINDER  MODEL  .  94 

REFERENCES . 96 

INITIAL  DISTRIBUTION  LIST  .  97 


IV 


ACKNOWLEDGMENTS 


I  express  my  sincere  thanks  and  appreciation  to  Dr. 
Young  S.  Shin  and  Dr.  Young  W.  Kwon  for  their  continued 
support,  guidance  and  encouragement  throughout  this 
research.  I  gratefully  acknowledge  the  continued  support  of 
the  underwater  shock  research  at  the  Naval  Postgraduate 
School  by  Doug  Bruder  and  Kent  Goering  of  the  Defense 
Nuclear  Agency. 

Many  others  not  mentioned  by  name  provided  welcome 
assistance.  To  you,  please  forgive  the  omission.  I  thank 
all  of  you. 

I  dedicate  this  work  to  my  wife  Marianne  and  my 
children,  Kyle  and  Hillary,  for  their  support  and 
understanding  over  the  duration  of  this  project. 


v 


I. 


INTRODUCTION 


The  use  of  USA/DYNA3D  for  analyzing  structures 
subjected  to  underwater  shock  has  been  shown  to  be  an 
effective  tool  in  predicting  the  response  of  these 
structures  and  the  resulting  damage  (Chisum  1992).  In  some 
cases  the  actual  deformation  of  a  cylinder  subjected  to  an 
underwater  shock  is  not  accurately  predicted.  One 
explanation  for  this  inaccuracy  is  that  actual  cylinders 
have  many  imperfections  (e.g.,  out  of  roundness,  thin 
sections,  voids,  etc.)  and  the  finite  element  modeling  of 
these  cylinders  often  do  not  take  into  account  the  initial 
imperfections  that  are  present.  By  introducing 
imperfections  in  the  position  of  the  node  points  generated 
by  a  finite  element  mesh  generator  it  is  possible  to  more 
accurately  model  the  geometry  of  an  actual  cylinder.  With  a 
more  accurate  model  of  the  actual  cylinder  the  finite 
element  numerical  analysis  of  the  resulting  damage  due  to 
underwater  shock  may  be  more  accurate  as  well. 
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II.  INITIAL  IMPERFECTIONS 


The  location  of  the  node  points  of  a  structure  modeled 
by  a  finite  element  mesh  generator  are  precisely  located  (to 
the  numerical  accuracy  of  the  modeling  code)  at  the 
positions  specified  in  the  inputs  to  the  mesh  generator 
algorithm.  Manufactured  structures  usually  have  many 
imperfections  that  are  not  accurately  modeled  by  many  finite 
element  mesh  generators.  Surveys  have  been  performed  to 
measure  the  imperfections  that  naturally  exist  in 
cylindrical  shell  structures  (Arbocz  1982).  These  data  have 
been  collected  into  data  banks,  and  Arbocz  suggests  that 
they  be  used  to  improve  design  criteria  for  buckling  of  thin 
shells.  The  imperfection  data  banks  show  that  imperfections 
have  characteristic  distributions  that  include  decreasing 
modal  amplitudes  with  increasing  mode  number.  Kirkpatrick 
(  1 9 8 9 )  found  that  by  introducing  initial  modal  imperf ect ions 
in  a  cylindrical  structure  subjected  to  blast  shock  loading 
that  the  resulting  numerical  analysis  agreed  much  closer  to 
the  experimental  data  than  the  analysis  without  these 
initial  imperfections. 

Initial  imperfections  were  introduced  into  the  location 
of  the  node  points  of  the  finite  element  model  using  a 


summation  of  modal  imperfections  expressed  as  the  cosine 
series  shown  in  Eq.  (1). 


N 

AR(  8)  =  £j4/>cos(«9  +  0/.) 

i»=2 


(1) 


where  AR  is  the  radial  imperfection,  0  is  the  angular 
position,  N  is  the  maximum  modal  contribution,  An  is  the 
modal  amplitude,  and  0*  is  a  random  modal  phase  shift.  The 

assumption  that  the  modal  phase  shift  is  a  random  variable 
is  reasonable  for  many  shells.  An  empirical  form  for  the 
modal  amplitude  is  expressed  as  shown  in  Eq.  (2). 


r 

n 


where  An  is  the  modal  amplitude  of  the  nth  modal 
imperfection,  n  is  the  mode  number  and  X  and  r  are 
coefficients  used  to  fit  the  available  data  for  shells  of  a 
given  construction.  The  modal  amplitude  can  be  modeled  as  a 
constant  percentage  of  the  shell  thickness.  For  many  of  the 
studies  in  this  investigation  the  modal  amplitude  is 
constant  as  shown  in  Eq.  (3)  at  1%  of  the  shell  thickness 
for  all  modes. 


/4n  =  0.01A  for  all  n 


(3) 
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Kirkpatrick  found  that  the  modal  amplitudes  in  Eq.  i4)  and 
(5)  accurately  modeled  the  imperfections  of  his  test 
cylinder . 

An  =  ().()5/i  n<  6 . (4) 

2  h 

A,  1  =  — r  n>l  . (  5  ) 

n' 

The  modal  imperfections  introduced  in  this  study  aie 
assumed  to  be  a  function  only  cf  the  circumfei  er.t  lal 
position  of  the  node  and  not  a  function  cf  the  axial 
position  of  the  node.  The  programs  used  tc  modify  the  node 
positions  are  shown  in  Appendix  A. 
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III.  TWO  DIMENSIONAL  INFINITE  CYLINDER  MODELS 


Fox  (1992)  developed  a  two  dimensional  infinite 
cylinder  model  to  validate  the  use  of  the  USA /DYNA3D  code 
for  use  in  underwater  shock  analysis.  This  model  was 
further  developed  into  the  full  symmetry  two  dimensional 
infinite  cylinder  model  shown  in  Figure  1.  This  model  has 
40  elements  around  the  circumf erence  of  the  cylinder,  radius 
of  6  inches,  shell  thickness  of  0.06  inches  and  an  element 
length  of  0.006  inches.  The  same  restriction  on  the  length 
of  the  elements  found  for  the  model  Fox  developed 
( length/ radius  =  0 . 001 )  was  also  true  for  this  model.  Fox 
stated  that  longer  elements  resulted  in  "residual  three 
dimensional  effects"  that  caused  oscillations  in  the 
numerical  results.  Symmetric  boundary  conditions  were 
imposed  on  both  ends  of  the  cylinder  to  model  an  infinite 
length  cylinder.  This  model  and  all  subsequent  models 
(unless  specifically  noted)  use  Belytschko-Tsay  (Belytschko 
1984)  shell  elements  for  their  numerical  efficiency  and  to 
eliminate  the  problems  of  using  a  solid  element  to  model  a 
thin  walled  structure.  For  details  of  the  model  input 
parameters  see  Appendix  B.  The  shell  material  is  modeled  as 
mild  steel,  an  elastic-plastic  material  with  a  Young’s 
modulus  of  2 . 9x1 07  psi,  a  Poisson's  ratio  of  0.3,  a 
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hardening  modulus  of  5100  psi  and  a  yield  stress  of  32,000 
psi.  For  details  of  the  material  models  used  in  this  study 
see  Appendix  C. 

A.  Perfect  Cylinder  Model 

This  cylinder  was  subjected  to  a  underwater  shock  wave 
modeled  as  a  plane  wave  with  an  amplitude  of  500  psi  and  a 
duration  of  1  millisecond  resulting  from  an  explosive  charge 
located  along  the  x-axis.  The  pressure  time  history  of  the 
shock  wave  is  shown  in  Figure  2.  The  early  damage  pattern 
of  the  cylinder  at  0.5  milliseconds  is  shown  in  Figure  3. 
The  front  face  of  the  cylinder  is  flattened  by  the  impacting 
plane  wave  approaching  from  the  right.  However,  Figure  4 
shows  that  the  final  damage  pattern  at  5  milliseconds  is 
very  different  than  this  early  deformation  pattern.  All 
plastic  deformation  has  been  completed  by  5  milliseconds. 
The  side  of  the  cylinder  facing  the  approaching  shock  wave 
shows  a  pronounced  protrusion.  This  pattern  was  unexpected 
based  on  experimental  results  which  have  usually  shown 
inward  deformations  of  a  cylindrical  shell  facing  the 
explosive  charge.  It  is  hypothesized  that  this  raised 
section  is  due  to  the  fact  that  the  modeled  cylinder  is  a 
perfect  cylinder  and  is  not  representative  of  a  true 
manufactured  cylinder.  The  manufactured  cylinder  will  have 
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4-Shock 

Wave 


0 . 5  ms  (deformations  scaled  up  100X) 

Figure  3.  Deformation  of  Two  Dimensional  Infinite  Cylinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1  ms 


4-Shock 

Wave 


5.0  ms  (deformations  scaled  up  10X) 

Figure  4.  Deformation  of  Two  Dimensional  Infinite  Cylinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1  ms 
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many  imperfections  in  the  shell  which  is  not  present  in  the 
computer  models. 

In  order  to  study  the  effect  of  mesh  size  on  the 
response  of  this  cylinder  another  model  was  developed  with 
64  elements  around  the  circumference  of  the  cylinder. 
Figure  5  shows  that  the  response  of  this  cylinder  is  similar 
to  that  for  the  40  element  cylinder.  There  is  still  a 
distinct  protrusion  of  the  shell  toward  the  direction  of  the 
approaching  shock  wave.  All  subsequent  cylinders  will  use 
40  elements  around  the  circumference  of  the  cylinder. 

B.  Imperfect  Cylinder  Models 

The  introduction  of  initial  imperfections  significantly 
changes  both  the  shape  and  magnitude  of  the  resulting 
deformation.  Introduction  of  a  6th  mode  imperfection,  as 
shown  in  Figure  6  (imperfection  magnitude  scaled  by  a  factor 
of  100)  with  a  modal  amplitude  of  5%  of  the  shell  thickness 
and  no  random  phase  shift,  results  in  the  damage  pattern 
shown  (deformations  scaled  by  factor  of  10)  in  Figure  7.  A 
mode  6  imperfection  was  chosen  because  of  the  ease  in 
identifying  its  distinctive  pattern.  Use  of  another  model 
imperfection  would  only  have  changed  the  shape  of  the  final 
deformation  not  any  conclusions  resulting  from  this  study. 
Comparison  of  the  final  damage  pattern  with  the  shape  of  the 
initial  imperfection  shows  that  the  damage  pattern  follows 
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<-Shock 

Wave 


5 . 0  ms  (deformations  scaled  up  10X) 

Figure  5.  Deformation  of  Two  Dimensional  Infinite  Cylinder 
Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1  ms 

Fine  Mesh 


Initial  Mode  6  Imperfection  (imperfections  scaled  up  100X) 

Figure  6.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1  ms,  Mode  6  Initial 
Imperfection  A6  =  0.05h  no  phase  shift 


4-Shock 

Wave 


5 . 0  ms  (deformations  scaled  up  10X) 

Figure  7.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1  ms,  Mode  6  Initial 
Imperfection  Ae  =  0.05h  no  phase  shift 
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the  initial  imperfection  resulting  in  a  final  damage  pattern 
looking  like  the  mode  6  initial  imperfection  magnified  due 
to  the  external  shock  wave  applied  to  the  cylinder.  The 
introduction  of  imperfections  caused  the  deformation  of  the 
cylindrical  shell  to  follow  the  initial  imperfection  pattern 
of  the  shell.  In  addition  to  changing  the  damage  pattern, 
the  initial  imperfection  also  causes  much  greater 
deflections  of  the  shell  for  the  same  shock  wave  intensity. 
Figure  8  shows  the  deformation  of  a  two  dimensional  cylinder 
with  an  initial  imperfection  5%  of  the  shell  thickness  with 
only  a  mode  10  contribution  and  no  random  phase  shift. 
Comparison  of  Figures  7  and  8  show  that  increasing  the  mode 
number  of  the  initial  imperfection  while  keeping  the 
amplitude  constant  at  5%  of  the  shell  thickness  results  in 
greater  damage  to  the  cylinder  for  the  same  shock  wave 
pressure.  This  is  due  to  the  greater  local  curvature  of  the 
shell  with  imperfections  of  higher  mode  number  than  for 
lower  mode  numbers.  Figure  9  shows  how  the  teal  internal 
strain  energy  (a  measure  of  the  damage  to  the  cylinder)  of 
the  cylinders  increases  as  the  mode  number  of  the  initial 
imperfection  increases  while  the  amplitude  of  the 
imperfection  is  constant  at  5%  of  the  shell  thickness  for 
the  same  shock  wave  pressure. 

Another  cylinder  was  modeled  with  imperfections 
containing  the  first  10  modal  imperfections,  as  shown  in 
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Internal  Strain  Energy  (in-lb{) 


5 . 0  ms  (deformations  scaled  up  3X) 

Figure  8.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1  ms,  Mode  10  Initial 
Imperfection  Aio  =  Q.05h  no  phase  shift 

2.0 

1.5 

1.0 


.5 


0 

0  .001  .002  .003  .004  .005 

Time  (sec) 


Figure  9.  Total  Internal  Strain  Energy  for  Cylinders 
subjected  to  a  Square  Pressure  Pulse  of  500  psi  for  1  ms 


Figure  10  with  modal  amplitudes  for  each  mode  of  1%  of  the 
shell  thickness  and  with  random  phase  shifts.  The  resulting 
damage  pattern  is  shown  in  Figure  11.  Again,  the  general 
shape  of  the  final  deformation  pattern  is  very  close  to  the 
shape  of  the  initial  imperfection  pattern. 

The  introduction  of  a  different  initial  imperfection, 
shown  in  Figure  12  with  initial  imperfections  magnified  by 
50  times,  using  the  modal  amplitude  of  equations  (4)  and 
(5),  results  in  the  final  deformation  pattern  shown  in 
Figure  13.  Again,  the  deformation  of  the  shell  of  the 
cylinder  follows  the  initial  imperfection  pattern. 
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Initial  Imperfection  (imperfections  scaled  up  200X) 


Figure  10.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1  ms,  Initial 
Imperfection  First  10  Modes  ^4 « =  0.01/j  with  Random  Phase  Shift 


<— Shock 
Wave 


5 . 0  ms  (deformations  scaled  up  2X) 

Figure  11.  Two  Dimensional  Infinite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  for  1  ms,  Initial 
Imperfection  First  10  Modes  A„  =  0.0lh  with  Random  Phase  Shift 
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IV.  THREE  DIMENSIONAL  RING  STIFFENED  INFINITE  CYLINDER 

MODELS 

The  three-dimensional  ring  stiffened  infinite  cylinder 
model  shown  in  Figure  14  was  developed  from  the  two- 
dimensional  infinite  cylinder  model.  This  model  has  a  total 
of  400  elements  with  40  elements  in  the  circumferential 
direction  and  10  elements  along  the  length.  Again  symmetric 
boundary  conditions  were  imposed  on  both  ends  of  the 
cylinder  to  model  an  infinite  cylinder.  This  chree- 
dimensional  model  has  stiffeners  0.12  inches  thick  and  1 
inch  deep  located  on  12  inch  spacing.  The  shell  is  0.06 
inches  thick.  For  details  in  the  model  input  parameters  see 
Appendix  D  The  shell  and  stiffeners  are  modeled  as  mild 
steel  with  the  same  properties  as  listed  earlier 

A.  Perfect  Cylinder  Model 

This  cylinder  was  subjected  to  the  same  plane  shock 
wave  as  the  previous  models.  The  early  deformation  of  the 
cylinder  at  0.2  milliseconds  (Figure  15)  shows  the  outer 
shell  pinching  in  on  either  side  of  the  stiffener  located  in 
the  center  of  the  cylinder.  Figure  16  shows  the  flattening 
of  the  face  of  the  cylinder  facing  the  explosive  charge  at 
0.2  milliseconds.  At  this  time  in  the  deformation  history 
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Figure  14. 


Three  Dimensional  Ring  Stiffened 
Cylinder  Model 
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Front  View  @  0.2  ms  (deformations  scaled  up  50X) 


Figure  15.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 


<— Shock 
Wave 


End  View  @  0.2  ms  (deformations  scale  up  50X) 

Figure  16.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 


there  is  no  evidence  of  the  formation  of  a  local  raised  area 
on  the  front  of  the  cylinder  facing  the  charge.  The  final 
deformation  pattern  shown  in  Figure  17  shows  pinching  of  the 
shell  on  either  side  of  the  stiffener.  This  pinching  effect 
was  also  noted  by  Chisum  (1992)  and  appears  to  be  present  in 
many  finite  element  simulations  of  underwater  shock  which 
model  test  cylinders  as  perfect  cylinders.  Again,  the  final 
deformation  pattern  of  the  cylinder  side  facing  the  charge 
has  a  local  protrusion  of  the  shell  material  toward  the 
explosive  charge  as  shown  in  Figure  18.  Again,  it  is 
hypothesized  that  this  damage  pattern  is  the  result  of  the 
model  cylinder  being  a  perfect  cylinder  while  test  cylinders 
will  always  have  imperfections  in  the  shell  of  the  cylinder. 

B.  Imperfect  Cylinder  Models 

The  introduction  of  initial  imperfections  significantly 
changes  both  the  shape  and  magnitude  of  the  final 
deformations  seen  in  this  cylinder  model.  The  5%  mode  6 
imperfection  shown  in  Figure  19  was  introduced  into  this 
model.  The  resulting  deformation  due  to  a  500  psi  1 
millisecond  plane  wave  pressure  pulse  is  shown  in  Figures  20 
through  23.  The  initial  imperfection  shape  is  clearly 
evident  in  the  final  damage  pattern.  However  the  shape  of 
the  cylinder  at  0.2  milliseconds  does  not  show  any  of  the 
initial  imperfection.  There  is  insufficient  time  elapsed 
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Front  View  @  5.0  ms 

Figure  17.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1  ms 


e-Shock 

Wave 


End  View  @  5.0  ms 

Figure  18.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1  ms 
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Mode  6  Initial  Imperfection  At,  =  0.05h  no  phase  shift 
(imperfections  scaled  up  100X) 

Figure  19.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  Subjected  to  a  Square  Pressure 

Pulse  of  500  psi  for  1  ms 


Front  View  @  0.2  ms  (deformations  scaled  up  50X) 


Figure  20.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Imperfection 
A6  =  0.05h  no  phase  shift,  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 


Shock 

Wave 


End  View  @  0.2  ms  (deformations  scaled  up  50X) 


Figure  21.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Imperfection 
A*  =  0.Q5h  no  phase  shift,  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 
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Front  View  @  5.0  ms 


Figure  22.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Imperfection 
A6  =  0.05/i  no  phase  shift,  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 


Shock 

Wave 


End  View  @  5.0  ms 


Figure  23.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Infinite  Length  Cylinder  with  a  Mode  6  Initial  Imperfection 
>46  =  0.05/?  no  phase  shift,  Subjected  to  a  Square  Pressure 
Pulse  of  500  psi  for  1  ms 


for  the  deformations  to  grow  large  enough  to  show  this 
initial  imperfection.  However  the  pinching  of  the  shell  on 


either  side  of 

the 

stiffener  has  already  begun  at 

0.2 

milliseconds . 

The 

final  deformation 

still  shows 

this 

pinching  effect 

however  the  magnitude 

of  this  type 

of 

deformation  is  reduced  by  the  initial  imperf ection. 

With  the  introduction  of  an  initial  imperfection  of  the 
first  10  modes  with  modal  amplitudes  of  1%  of  the  shell 
thickness  and  with  random  phase  shifts  as  shown  in  Figure  24 
the  resulting  damage  pattern  again  shows  that  the 
deformation  of  the  shell  will  preferentially  fellow  the 
initial  imperfections.  The  resulting  damage  pattern  is  very 
different  from  the  perfect  case  (Figures  25  and  26).  Most 
noteworthy  is  the  elimination  of  the  pinch  in  the  outer 
shell  on  either  side  of  the  stiffener.  In  addition  the 
magnitude  of  the  deformation  in  the  outer  shell  is  much 
greater  than  for  the  perfect  cylinder  case.  The 
introduction  of  these  imperfections  results  in  sites  where 
the  cylinder  preferentially  deforms  during  the  shock 
pressure.  Thus  when  the  cylinder  deforms  the  deformation 
follows  the  initial  imperfections  resulting  in  the  final 
shape  of  the  cylinder  looking  like  the  initial  imperfection 
shape . 
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Infinite  Length 
of  500  psi  for 


Front  View  @  5.0  ms 


Figure  25.  Three  Dimensional  Ring  Stiffened  Infinite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1  ms,  Initial  Imperfection  First  10  Modes  /U  =  0.01/i  with 

Random  Phase  Shift 


<— Shock 
Wave 


End  View  @  5 . . 0  ms 


Figure  26.  Three  Dimensional  Ring  Stiffened  Infinite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1  ms,  Initial  Imperfection  First  10  Modes  An  =  0.0\h  with 

Random  Phase  Shift 
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C.  Parametric  Study  on  the  Effect  of  Cylinder  Geometry 


A  parametric  study  changing  basic  cylinder  geometry  was 

done  to  see  the  effect  of  these  changes  on  the  response  of  a 

three  dimensional  ring  stiffened  infinite  cylinder  to 

underwater  shock.  The  parameters  were  varied  as  follows: 

Shell  Thickness  0.06,  0.12  and  0.24  inches 
Stiffener  Spacing  12.0,  8.0  and  6.0  inches 
Diameter  8.0,  12.0  and  18.0  inches 

This  geometric  variation  covers  the  range  of  shell  thickness 

to  radius  ratios  from  0.0067  to  0.06  and  stiffener  spacing 

to  radius  ratios  from  0.67  to  3.0. 

Increasing  the  shell  thickness  strengthened  the 
cylinders  and  in  order  to  obtain  significant  permanent 
deformation  of  the  shell  it  was  necessary  to  increase  the 
amplitude  of  the  input  shock  wave.  ^he  amplitude  of  the 
shock  wave  was  adjusted  to  cause  damage  to  the  cylinder 
without  completely  collapsing  the  cylinder.  This  input 
shock  wave  pressure  ranged  from  500  psi  for  1  ms  for  the 
thinnest  shell  to  2400  psi  for  1  ms  for  the  thickest  shell. 

The  resulting  damage  patterns  for  a  perfect  cylinder 
8.0  inches  in  diameter  (shell  thickness  to  radius  ratios 
from  0.015  to  0.06  and  stiffener  spacing  to  radius  ratios 
from  1.5  to  3.0)  are  shown  in  Figures  27  and  28.  The  shock 
wave  pressure  for  the  0.06  inch  shell  thickness  is  600  psi 
for  1  ms,  for  the  0.12  inch  shell  thickness  the  pressure  is 
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1200  psi  for  1  ms  and  for  the  0.24  inch  shell  thickness  the 
pressure  is  2400  psi  fcr  1  ms. 

The  front  view  shows  that  increasing  the  shell 
thickness  causes  the  pinching  of  the  shell  near  the 
stiffener  to  decrease  but  it  is  not  entirely  eliminated  at  a 
shell  thickness  of  0.24  inches.  Decreasing  the  stiffener 
spacing  to  8.0  inches  from  12.0  inches  also  -v\  little 
effect  on  the  pinching  of  the  shell  but  at  a  stiffener 
spacing  of  6.0  inches  there  is  very  little  pinching  of  the 
shell  near  the  stiffener. 

The  end  view  shows  that  at  a  shell  thickness  of  0.06 
inches  there  is  some  damage  to  the  front  and  back  of  the 
cylinder  but  most  of  the  damage  occurs  to  the  top  and  bottom 
of  the  cylinder.  For  a  stiffener  spacing  of  12.0  inches  and 
a  shell  thickness  of  0.12  inches  there  is  still  significant 
damage  to  the  cylinder  at  the  top  and  bottom  but  there  is 
additional  damage  to  the  front  and  back  of  the  cylinder. 
Finally  at  a  stiffener  spacing  of  12.0  inches  and  a  shell 
thickness  of  0.24  inches  most  of  the  damage  to  the  cylinder 
is  at  the  back  with  some  damage  to  the  front  and  relatively 
little  damage  to  the  top  and  bottom.  Reducing  the  stiffener 
spacing  also  effects  the  damage  pattern  by  causing  more 
damage  to  the  front  of  the  cylinder  and  less  to  the  top  and 
bottom  of  the  cylinder  than  for  the  12.0  inch  stiffener 
spacing . 
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The  resulting  damage  patterns  for  a  perfect  cylinder 
12.0  inches  in  diameter  (shell  thickness  to  radius  ratios 
from  0.01  to  0.04  and  stiffener  spacing  to  radius  ratios 
from  1.0  to  2.0)  are  shown  in  Figures  29  and  30.  The  shock, 
pressure  for  the  0.06  inch  shell  thickness  is  500  psi  for  1 
ms,  for  the  0.12  inch  shell  thickness  the  pressure  is  1000 
psi  for  1  ms  and  for  the  0.24  inch  shell  thickness  the 
pressure  is  1900  psi  for  1  ms. 

At  a  shell  thickness  of  0.06  inches  the  pinching  of  the 
shell  occurs  for  stiffener  spacing  of  12.0.  8.0  and  6.0 

inches.  At  a  stiffener  spacing  of  12.0  inches  increasing 
the  shell  thickness  causes  the  pinching  of  the  shell  to 
decrease  until  it  is  eliminated  for  the  shell  thickness  of 
0.24  inches.  At  stiffener  spacing  of  less  than  12.0  inches 
increasing  the  shell  thickness  causes  the  pinching  of  the 
shell  to  disappear  for  both  the  0.12  and  0.24  inch  shell 
thickness . 

The  end  view  shown  that,  for  shell  thickness  of  0.06 
and  0.12  inches,  at  a  diameter  of  12.0  inches  there  is  a 
protrusion  of  the  shell  toward  the  shock  wave  approach 
direction.  Increasing  the  shell  thickness  to  0.24  inches, 
for  the  three  stiffener  spacing  investigated,  causes  the 
protrusion  to  be  eliminated.  Decreasing  the  stiffener 
spacing  causes  a  decrease  in  damage  to  the  cylinder  and  some 
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minor  shifting  of  the  damage  pattern  but  the  basic  shape  of 
the  damages  pattern  remains  constant. 

The  resulting  damage  for  a  perfect  cylinder  18.0  inches 
in  diameter  (shell  thickness  to  radius  ratios  from  0.0067  to 
0.0267  and  stiffener  spacing  to  radius  ratios  from  0.67  to 
1.33)  is  shown  in  Figures  31  and  32.  The  shock  pressure  for 
the  0.06  inch  shell  thickness  is  500  psi  for  1  ms,  for  the 
0.12  inch  shell  thickness  the  pressure  is  1000  psi  for  1  ms 
and  for  the  0.24  inch  shell  thickness  the  pressure  is  1700 
psi  for  1  ms. 

In  all  cases,  for  the  18.0  inch  diameter  perfect 
cylinder  there  was  no  pinching  of  the  shell  near  the 
stiffener.  This  cylinder  is  structurally  much  less  stiff 
than  the  smaller  diameter  cylinders  thus  allowing  the  shell 
to  collapse  more  readily  without  pinching  near  the 
stiffener . 

From  the  end  view  it  can  be  seen  that  the  protrusion  of 
the  shell  toward  the  explosion  source  occurs  in  all  cases. 
For  the  thicker  shell  this  protrusion  is  not  as  sharp  as  for 
the  thinner  shell.  Decreasing  the  stiffener  spacing  does 
not  greatly  affect  the  damage  patterns  but  does  cause  the 
damage  levels  to  decrease. 

The  damage  pattern  for  a  12.0  inch  diameter  cylinder, 
with  a  mode  6  initial  imperfection,  an  amplitude  of  5%  of 
the  shell  thickness  and  no  phase  shift,  is  shown  in  Figures 
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33  and  34.  The  shock  wave  pressure  varies  from  LOO  psi  for 
1  ms  for  the  0.06  inch  shell  thickness  to  1000  psi  for  1  ms 
for  the  0.12  inch  shell  thickness  and  finally  to  1900  psi 
for  1  ms  for  the  0.24  inch  shell  thickness.  These  shock 
wave  pressure  are  the  same  as  for  the  perfect  cylinder  cases 
at  a  cylinder  diameter  of  12.0  inches. 

The  front  view  shows  the  pinching  of  the  shell  near  the 
stiffener  only  occurs  for  the  0.06  inch  shell  thickness. 
For  thicker  shells  the  pinching  is  eliminated.  The  end 
views  show  that  the  initial  mode  6  imperfection  is  seen  in 
the  final  deformation  pattern  of  the  shell  of  the  cylinder 
in  all  cases.  Changing  the  shell  thickness  and  the 
stiffener  spacing  has  very  little  effect  on  the  damage 
pattern.  The  damage  level  is  reduced  as  the  stiffener 
spacing  is  reduced  and  this  is  expected  as  the  cylinder 
becomes  stiffer  with  decreasing  stiffener  spacing. 

The  damage  to  an  8.0  inch  cylinder,  with  initial 
imperfections  consisting  of  the  first  10  mode  shapes  with  an 
amplitude  of  1%  of  the  shell  thickness  and  witn  random  phase 
shifts  of  the  modes  is  shown  in  Figures  35  and  36.  In  these 
cases,  the  shock  wave  pressures  had  to  be  reduced  from  the 
pressures  used  in  the  perfect  cylinders  at  the  same 
diameter,  because  the  cylinders  completely  collapsed  at  the 
pressures  used  in  the  perfect  cylinder  cases.  The  shock 
wave  pressures  were  500  psi  for  1  ms  for  the  0.06  inch  shell 
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thickness,  1000  psi  for  1  ms  for  the  0.12  inch  shell 
thickness  and  2250  psi  for  1  ms  for  the  0.24  inch  shell 
thickness.  At  this  diameter  the  cylinder  was  weakened  to  a 
greater  extent,  than  for  the  larger  diameter  cases,  by  the 
initial  imperfections  introduced  into  the  shell  geometry. 

The  front  view  of  the  cylinder  shows  that  there  is  very 
little  pinching  of  the  shell  of  the  cylinder  near  the 
stiffeners.  The  end  view  shows  the  damage  is  much  more 
localized  to  a  smaller  portion  of  the  shell  than  for  the 
larger  diameter  cylinders.  The  damage  pattern  is  fairly 
consistent  with  only  changes  in  the  magnitude  of  the  damage 
among  the  cases  studied.  In  addition  the  initial 
imperfection  pattern  is  still  seen  in  the  final  damage 
pattern  but  the  match  is  much  less  pronounced  than  for  the 
later  cases  studied. 

For  a  cylinder  diameter  of  12.0  inches  with  initial 
imperfections  consisting  of  the  first  10  modes  shapes  with 
an  amplitude  of  1%  of  the  shell  thickness  with  random  phase 
shifts  of  the  modes,  the  damage  patterns  are  much  more 
consistent.  Figures  37  and  38  shown  the  resulting  damage 
patterns.  At  this  diameter  the  shock  wave  pressure  is  the 
same  as  for  the  perfect  cylinder  cases  at  the  same  diameter 
with  pressures  of  500  psi  for  1  ms  for  the  0.06  inch  shell 
thickness,  1000  psi  for  1  ms  for  the  0.12  inch  shell 
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thickness  and  1900  psi  for  1  ms  for  the  0.24  inch  shell 
thickness . 

The  front  view  shows  a  similar  damage  pattern  for  all 
cases.  The  shell  very  distinctly  dishes  between  the 
stiffeners  in  a  manner  similar  to  that  seen  in  many 
underwater  shock  test  cylinders.  Changing  the  shell 
thickness  and  the  stiffener  spacing  has  very  little  effect 
on  the  damage  pattern  as  seen  from  the  front  of  the 
cylinder.  The  end  view  shows  again  that  in  all  cases  shown 
the  damage  pattern  remains  very  consistent  and  only  the 
damage  levels  change  with  changes  in  the  stiffener  spacing. 
The  0.24  inch  shell  thickness  does  show  some  shifting  of  the 
damage  pattern  but  the  same  basic  pattern  is  still  evident. 

Finally,  Figures  39  and  40  show  the  deformation  of  an 
18.0  inch  diameter  cylinder  with  initial  imperfections 
consisting  of  the  first  10  mode  shapes,  each  with  an 
amplitude  of  1%  of  the  shell  thickness  with  random  phase 
shift.  Again,  the  shock  wave  pressure  was  the  same  as  for 
the  perfect  cylinder  cases  with  a  pressure  of  500  psi  for  1 
ms  for  the  0.06  inch  shell  thickness,  1000  psi  for  1  ms  for 
the  0.12  inch  shell  thickness  and  1700  psi  for  1  ms  for  the 
0.24  inch  shell  thickness. 

From  the  front  view  it  can  be  seen  that  the  shell  of 
the  cylinder  has  deformed  by  dishing  between  the  stiffeners 
in  all  cases.  Decreasing  the  stiffener  spacing  reduced  the 
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magnitude  of  the  deformation  but  did  not  change  the  shape  of 
the  deformation.  Increasing  the  shell  thickness  had  little 
effect  on  the  deformation  pattern.  The  end  views  of  the 
cylinder  all  show  the  same  basic  deformation  pattern. 
Increasing  the  shell  thickness  had  some  effect  on  the  damage 
pattern  but  the  same  deformation  pattern  is  present  in  all 
the  deformed  cylinders. 
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V.  THREE  DIMENSIONAL  RING  STIFFENED  FINITE  LENGTH  CYLINDER 

From  the  three-dimensional  infinite  cylinder  a  three 
dimensional  finite  length  cylinder  was  developed  as  shown  in 
Figure  41.  This  cylinder  is  three  feet  long,  1  foot  in 
diameter  with  two  stiffeners  evenly  spaced  12  inches  apart. 
The  shell  of  the  cylinder  is  mild  steel  0.06  inches  thick. 
The  stiffeners  are  mild  steel  0.12  inches  thick  and  1  inch 
deep.  The  endplate  is  HY-100  steel  0.25  inches  thick.  The 
HY-100  steel  is  modeled  as  an  elastic-plastic  material  with 
a  Young's  modulus  of  2.9xl07  psi,  a  Poisson's  ratio  of  0.3, 
a  hardening  modulus  of  5020  psi  and  a  yield  stress  of 
108,000  psi.  This  cylinder  is  modeled  as  a  half  cylinder 
with  a  plane  of  symmetry  perpendicular  to  the  axis  of 
rotation.  The  model  has  40  elements  in  the  circumferential 
direction  and  15  elements  in  the  axial  direction  for  a  total 
of  600  elements  and  921  nodes.  The  use  of  a  half  symmetry 
model  with  a  symmetric  boundary  condition  results  in  a 
smaller  number  of  elements  and  greater  computational 
efficiency  for  the  finite  element  analysis.  For  details  of 
the  model  input  parameters  see  Appendix  E. 
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Figure  41.  Three  Dimensional  Ring  Stiffened  Finite  Length 

Cylinder  Model 
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A.  Perfect  Cylinder  Models 

This  cylinder  was  subjected  to  the  same  planar  shock 
wave  as  in  the  previous  models.  The  resulting  damage 
pattern  for  a  perfect  cylinder  can  be  seen  in  Figures  42 
through  49.  Again  the  side  of  the  cylinder  facing  the 
planar  shock  wave  shown  a  local  raised  area  that  was  also 
seen  in  the  previous  models.  In  addition  there  is  a  pinch 
in  the  shell  near  the  endplates  and  on  either  side  of  the 
stiffeners.  This  pinch  was  also  noted  by  Chisum  (1992) 
during  his  analysis  of  DNA  models  used  for  underwater  shock 
testing.  Most  of  the  deformation  of  the  cylinder  shell 
occurs  between  1.0  and  2.0  milliseconds  (Figures  44  through 
47).  During  the  time  from  1.0  to  2.0  milliseconds  this 
kinetic  energy  is  transferred  to  strain  energy  causing 
deformation  of  the  shell.  after  2.0  milliseconds  the  most 
of  the  plastic  deformation  has  occurred  and  the  cylinder 
shell  then  deforms  only  through  elastic  vibration  (Figures 
48  and  49)  .  Figure  50  is  a  plot  of  the  kinetic  energy  of 
the  cylinder  shell  as  a  function  of  time.  The  shell  has  its 
maximum  kinetic  energy  at  1  ms  and  then  this  energy  rapidly 
decreases  as  it  is  dissipated  into  the  surrounding  medium 
and  into  the  shell  material  as  strain  energy.  Figure  51 
shows  the  strain  history  of  the  element  on  the  front  of  the 
cylinder  exactly  in  the  middle  of  the  stiffeners.  From  this 
plot,  at  1  millisecond  the  hoop  strain  is  -5590  microstrain. 
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Front  View  @  0.2  ms  (deformations  scaled  up  50X) 


Figure  42.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  0.2  ms  (deformations  scaled  up  50X) 

Figure  43.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 


Front  View  @  1.0  ms  (deformations  scaled  up  2X) 

Figure  44.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  1.0  ms  (deformations  scaled  up  2X) 

Figure  45.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  2.0  ms  (deformations  scaled  up  2X) 

Figure  46.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  2.0  ms  (deformations  scaled  up  2X) 

Figure  47.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  5.0  ms  (deformations  scaled  up  2X) 

Figure  48.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  (deformations  scaled  up  2X) 

Figure  49.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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at  2  milliseconds  the  hoop  strain  is  -15000  microstrain  and 
at  5  milliseconds  the  hoop  strain  is  -14600  microstrain. 
The  strain  at  5  milliseconds  is  permanent  plastic  strain. 
Most  of  this  plastic  strain  occurred  between  1  and  2 

milliseconds  with  very  little  additional  strain  occurring 

after  2  milliseconds.  This  confirms  the  previously  noted 

timing  of  the  deformation  of  the  shell  of  the  cylinder. 
Most  of  the  damage  to  the  shell  occurs  between  1  and  2 

milliseconds . 

B.  Imperfect  Cylinder  Models 

The  addition  of  initial  imperfections  greatly  changes 
the  resulting  deformation  pattern  of  the  shell  of  the 
cylinder.  The  addition  of  a  mode  6  imperfection  (shown  in 
Tigure  52)  with  a  modal  amplitude  of  5%  of  the  shell 
thickness  and  with  no  random  phase  shift  results  in  the 
final  damage  pattern  shown  in  Figure  53  and  54.  The  end 
view  clearly  shows  the  strong  effect  of  the  mode  6  initial 
imperfection  causing  the  final  deformation  to  follow  this 
initial  imperfection.  The  front  view  of  the  cylinder  shows 
that  the  pinching  of  the  shell  near  the  endplates  and 
stiffeners  is  very  much  reduced  due  to  the  introduction  of 
this  imperfection. 

An  investigation  into  the  effect  of  the  modal  amplitude 
on  the  final  deformation  pattern  showed  that  very  small 
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Front  View  @  5.0  ms 

Figure  53.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1  ms,  Mode  6  Initial  Imperfection  /46  =  0.05/i  no  phase  shift 
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End  View  @  5.0  ms 

Figure  54.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Square  Pressure  Pulse  of  500  psi  for 
1  ms.  Mode  6  Initial  Imperfection  >U  =  0.05/i  no  phase  shift 
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modal  amplitudes  have  a  significant  effect  on  the  response 
of  the  cylinder.  An  initial  imperfection  consisting  of  the 
first  10  mode  shapes  with  varying  modal  amplitudes  and 
random  phase  shifts  is  shown  in  Figure  55.  The  shell 
deformation  again  follows  the  initial  imperfections  present 
in  the  shell  structure.  As  the  modal  amplitudes  increase 
the  deformation  of  the  shell  of  the  cylinder  also  increases. 
At  a  modal  amplitude  of  0.01%  of  the  shell  thickness 
(Figures  56  and  57)  the  front  of  the  cylinder  is  distinctly 
raised  toward  the  explosive  charge  as  in  the  perfect 
cylinder.  But,  it  can  also  be  seen  from  the  end  view  that 
the  deformation  pattern  even  at  this  magnitude  of 
imperfection  has  changed  from  the  pattern  for  the  perfect 
cylinder.  At  a  model  amplitude  of  0.05%  (Figures  58  and  59) 
this  alteration  of  the  pattern  is  quite  obvious.  The 
pinching  of  the  shell  of  the  cylinder  persists  until  higher 
modal  amplitudes  are  reached  (Figures  60  through  63). 
Finally  at  a  modal  amplitude  of  0.5%  of  the  shell  thickness 
(Figures  62  and  63)  the  pinching  of  the  shell  is  eliminated. 
Increasing  the  amplitude  further  results  in  increased  damage 
to  the  shell  due  to  the  shock  wave  (Figures  64  through  67). 
It  should  be  noted  that  an  amplitude  of  0.5%  of  the  shell 
thickness  represents  a  modal  amplitude  of  0.0003  inches  for 
this  model.  This  magnitude  of  imperfections  would  almost 
certainly  be  present  in  a  test  cylinder. 
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Front  View  @  5.0  ms  A*  =  0.000 Ih  with  Random  Phase  Shift 

Figure  56.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  >4<.  =  0.0001/j  with  Random  Phase  Shift 

Figure  57.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  5.0  ins  A*  =  0.0005h  with  Random  Phase  Shift 

Figure  58.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  /4«  =  0.0005A  with  Random  Phase  Shift 

Figure  59.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  5.0  ms  An  =  0.001/?  with  Random  Phase  Shift 

Figure  60.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  A«  =  0.00I/i  with  Random  Phase  Shift 

Figure  61.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  5.0  ms  An  =  0.005/?  with  Random  Phase  Shift 

Figure  62.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  An  =  0.QQ5h  with  Random  Phase  Shift 

Figure  63.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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Front  View  @  5.0  ms  A*  =  0.01/?  with  Random  Phase  Shift 

Figure  64.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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End  View  @  5.0  ms  An  =  0.01A  with  Random  Phase  Shift 

Figure  65.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 


68 


Front  View  @  5.0  ms  A*  =  0.Q5h  with  Random  Phase  Shift 

Figure  66.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 


Shock— » 
Wave 


End  View  @  5.0  ms  A»  =  0.05/i  with  Random  Phase  Shift 

Figure  67.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 

of  500  psi  for  1  ms 
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The  development  of  a  raised  section  on  the  cylinder 
facing  the  shock  wave  was  unexpected.  With  the  previous 
perfect  cylinder  the  shock  wave  impacted  the  cylinder  along 
the  line  of  nodes  at  the  front  of  the  cylinder.  The 
cylinder  was  then  rotated  so  that  the  shock  wave  would 
impact  the  cylinder  along  a  line  of  elements.  The  resulting 
deformation  pattern  is  shown  in  Figure  68.  There  is  still  a 
local  raised  section  on  the  cylinder  facing  the  shock  wave 
direction.  Thus,  it  is  felt  that  the  shot  geometry  of  the 
finite  element  mesh  is  not  a  contributing  factor  to  the 
development  of  this  raised  section. 

In  addition,  a  very  specific  imperfection  was 
introduced  in  the  five  nodes  facing  the  explosive  charge 
centrally  located  in  the  midbay  of  the  cylinder.  This 
imperfection  was  an  inward  imperfection  of  the  node 
locations  with  an  amplitude  of  1%  of  the  shell  thickness 
(0.006  inches).  The  damage  pattern  resulting  from  a  500  psi 
1  millisecond  shock  wave  is  shown  in  Figure  69. 

The  use  of  Hughes-Liu  shell'  element  (Hughes  1981)  vice 
the  Belytschko-Tsay  shell  element  was  investigated.  The 
Hughes-Liu  shell  element  does  not  use  many  of  the 
simplifying  assumptions  used  in  formulating  the  Belytschko- 
Tsay  shell  element.  As  a  result,  the  use  of  the  Hughes-Liu 
shell  element  requires  longer  computational  times  that  may 
not  be  acceptable  in  certain  models.  Figures  70  and  71  show 
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End  View  @  5.0  ms 

Figure  68.  Deformation  of  Three  Dimensional  Ring  Stiffened 
Finite  Length  Cylinder  Subjected  to  a  Square  Pressure  Pulse 
of  500  psi  for  1  ms.  Cylinder  Rotated  4.5  Degrees 
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End  View  @  5.0  ms 

Figure  69.  Three  Dimensional  Finite  Cylinder  Subjected  to  a 
Square  Pressure  Pulse  of  500  psi  1  ms,  Initial  Imperfection 
Front  Midbay  Nodes  1%  Offset  Radially  Inward 
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the  resulting  deformation  to  a  model  using  Hughes-Liu  shell 
elements.  This  model  is  identical  to  the  previous  models 
with  the  exception  of  the  use  of  the  Hughes-Liu  shell 
element.  This  model  has  initial  imperfections  with  the 
first  1C  modes  and  modal  amplitudes  of  1%  of  the  shell 
thickness  and  random  phase  shifts.  Comparison  of  Figures  70 
and  71  with  Figures  64  and  65  shows  that  the  Hughes-Liu  and 
the  Belytschko-Tsay  shell  elements  have  the  same  final 
deformations.  In  this  case  the  use  of  the  numerically  more 
complicated  Hughes-Liu  shell  element  is  not  warranted. 

C.  Exponential  Decay  Shock  Wave 

The  use  of  an  explosive  shock  wave  modeled  as  a  square 
wave  is  representative  of  a  very  large  explosion  at  a  great 
standoff  distance  from  the  cylinder.  In  order  to  model  a 
smaller  explosion  close  to  the  cylinder  a  pressure  profile 
corresponding  to  that  produced  by  40  pounds  of  PETN  at  a 
standoff  distance  of  30  feet  was  used.  The  pressure  time 
history  of  the  resulting  shock  wave  is  shown  in  Figure  72. 
This  explosion  produces  a  peak  pressure  of  1839  psi  and  an 
exponentially  decaying  pressure  history  with  a  time  constant 
of  0.311  milliseconds.  The  resulting  damage  to  the  cylinder 
from  this  shock  wave  is  shown  in  Figure  73  through  75.  From 
the  offset  view  it  can  be  seen  that  the  response  of  the 
cylinder  is  different  from  that  to  the  previous  plane  wave. 
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Front  View  @  5.0  ms 


Figure  70.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Hughes-Liu  Shell  Elements  Subjected  to  a  Square 
Pressure  Pulse  of  500  psi  for  1  ms,  Initial  Imperfection 
First  10  Modes  A*  =  0.0l/i  with  Random  Phase  Shift 
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End  View  @  5.0  ms 

Figure  71.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Hughes-Liu  Shell  Elements  Subjected  to  a  Square 
Pressure  Pulse  of  500  psi  for  1  ms.  Initial  Imperfection 
First  10  Modes  An  =  0.01/1  with  Random  Phase  Shift 
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Figure  72.  Pressure  Time  History  of  Shock  Wave  from  40  its 
PETN  with  a  30  foot  Standoff  and  1  ms  Surface  Cutoff  Time 


Front  View  @  5.0  ms 

Figure  73.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 

Foot  Standoff 
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The  small  standoff  distance  causes  the  shock  wave  to  be  a 
spherical  wave  at  the  cylinder.  The  midbay  still  has  a 
protrusion  of  the  shell  material  toward  the  explosive  charge 
however  the  endbay  shell  material  has  a  depression  of  the 
shell  facing  the  charge.  Also  evident  from  the  front  view 
is  the  severe  pinching  of  the  shell  near  the  endplate  and  on 
either  side  of  the  stiffener.  This  pinching  was  very  severe 
near  the  endplate. 

The  introduction  of  a  5%  mode  6  imperfection  to  the 
model  (Figure  76)  results  in  the  deformations  shown  in 
Figure  77  and  78.  Again,  it  is  evident  that  the  damage 
pattern  of  the  cylinder  due  to  the  explosive  shock  clearly 
followed  the  initial  imperfection.  The  pinching  of  the 
shell  material  near  the  endplate  and  stiffeners  was  reduced 
but  not  eliminated. 

Introducing  an  imperfection  of  the  first  10  mode  shapes 
with  a  modal  amplitude  of  1%  of  the  shell  thickness  with 
random  phase  shifts  (Figure  79)  results  in  the  final  damage 
pattern  shown  in  Figures  80  and  81.  The  deformation  of  the 
shell  of  the  cylinder  was  again  influenced  by  the  initial 
imperfection  pattern.  The  local  pinching  of  the  shell 
material  on  either  side  of  the  stiffener  was  eliminated  and 
the  pinching  near  the  endplate  was  reduced  in  magnitude. 
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Initial  Imperfection  (imperfections  scaled  up  100X) 

Figure  76.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Mode  6  Imperfection  At  =  0.05h  no  phase 

shift 
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Front  View  @  5.0  ms 

Figure  77.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Mode  6  Imperfection  >46  =  0.05/?  no  phase 

shift 
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End  View  @  5.0  ms 

Figure  78.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Mode  6  Imperfection  >46  =  0.05/?  no  phase 

shift 
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Initial  Imperfection  (imperfections  scaled  up  200X) 

Figure  79.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Imperfection  First  10  modes  A « =  0. 0 l/i 

with  Random  Phase  Shift 


Front  View  @  5.0  ms 

Figure  80.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Imperfection  First  10  modes  /l*  =  0.01/? 

with  Random  Phase  Shift 
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End  View  @  5.0  ms 

Figure  81.  Three  Dimensional  Ring  Stiffened  Finite  Length 
Cylinder  Subjected  to  a  Pressure  Pulse  from  40  lbs  PETN  30 
Foot  Standoff,  Initial  Imperfection  First  10  modes  An  =  0.01/? 

with  Random  Phase  Shift 
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VI.  SUMMARY  AND  CONCLUSIONS 


The  response  of  test  cylinders  to  underwater  explosive 
shock  is  a  complicated  function  of  many  factors.  One 
important  factor  that  is  often  overlooked  in  the  modeling  of 
these  underwater  shock  phenomena  is  the  initial 
imperfections  that  are  present  in  the  test  cylinders.  The 
response  of  model  cylinders  subjected  to  simulated 
underwater  explosive  shock  has  been  shown  to  be  very 
dependent  to  the  initial  imperfections  introduced  into  the 
cylinder  models.  The  introduction  of  these  imperfections 
not  only  caused  the  shape  of  the  shell  of  the  cylinder  (as 
viewed  from  the  end)  to  follow  the  shape  of  the  initial 
imperfection  but  it  has  also  changed  the  response  of  the 
shell  near  the  endplates  and  stiffeners  of  the  model 
cylinder.  The  pinching  of  the  shell  near  these  stiffeners 
and  endplates  was  greatly  reduced  or  eliminated  by  the 
introduction  of  initial  imperfections  in  the  model 
cylinders.  This  resulting  response  is  closer  to  the 
response  observed  in  test  cylinders  subjected  to  actual 
underwater  shock  loading.  In  addition,  the  initial 
imperfections  introduced  increased  the  magnitude  of  the 
shell  deformations  compared  to  the  deformations  for 
cylinders  modeled  as  perfect  cylinders. 
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Changes  in  the  cylinder  geometry  resulted  in  some 
changes  in  the  response  of  the  cylinder  to  underwater  shock. 
However  in  all  cases  the  introduction  of  initial 
imperfections  greatly  affected  the  response  of  the  cylinder 
as  compared  to  the  response  of  a  perfect  cylinder. 

If  the  initial  imperfections  of  a  test  cylinder  are 
known  and  introduced  into  a  model  cylinder  used  in  finite 
element  analysis,  the  results  of  the  finite  element  analysis 
may  more  closely,  simulate  the  actual  response  of  the  test 
cylinder . 
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APPENDIX  A 


FORTRAN  PROGRAMS  FOR  MODIFYING  INGRIDO  FILE  FOR  MODAL 

IMPERFECTIONS 


program  imp 

************************************************************** 

*  This  program  written  by  Donald  T.  Hooker  II  on  3/15/93  * 

*  Revised  5/12/93  * 

*  **•****★***•******•*•******************•*************•*,«■**•*•****** 

*  This  program  modifies  the  nodes  of  a  file  called  * 

*  " ingrido . raw" ,  which  is  an  edited  version  of  ingrido  in  * 

*  which  the  first  line  is  the  beginning  of  the  node  list.  * 

*  This  is  one  of  a  series  of  2  programs,  which  need  to  be  * 

*  run  in  order,  i.e.  "imp"  then  "convert".  These  programs* 

*  modify  the  node  positions  for  modal  imperfections  WITH  * 

*  OR  WITHOUT  random  amplitudes  and/or  phase  shifts.  The  * 

*  final  output  file  name  is  imp.dat.  * 


★ 


★ 


★ 


★ 

★ 

★ 

★ 

★ 


★ 


VARIABLE  DECLARATION 

a  -  angle  from  reference  axis  to  radial  position  of  * 
node  * 

al  thru  alO  -  weighting  factor  for  each  modal  * 

imperfection  * 

bl  -  ingrido  constants  * 

b2  -  ingrido  constants  * 

dr  -  change  in  radial  position  due  to  modal  * 

imperfections  * 

drl  thru  drlO  -  radial  position  change  for  each  modal  * 
imperfection  * 

h  -  thickness  of  shell  * 

n  -  node  number  * 

m  -  number  of  nodes  * 

pi  -  3.14159265359 

psl  thru  pslO  -  phase  shift  for  each  modal  * 

imperfection  * 

r  -  radial  position  of  node  * 

randamp  -  if  random  amplitude  is  desired  set  randamp=l* 
randphs  -  if  random  phase  is  desired  set  randphs=l  * 
x  -  x  coordinate  of  node  * 

y  -  y  coordinate  of  node  * 

z  -  z  coordinate  of  node  * 


************************************************************ 


double  precision  x,y , z , a , r , dr , h, pi 

double  precision  al , a2 , a3 , a4 , a5 , a6 , a7 , a8 , a9 , alO 

double  precision  drl , dr2 , dr3 , dr4 , dr5 , dr6 , dr7 , dr8 , dr9 , 
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&  drlO 

double  precision  psl , ps2 , ps3 , ps4 , ps5 , ps6 , ps7 , ps8 , ps9 , 

&  pslO 

integer  randamp, randphs 

************************************************************ 

*  OPEN  INPUT  AND  OUTPUT  FILES  * 

************************************************************ 

open ( 14 , f ile= 1 ingrido . raw ' ) 
open ( 15, f ile= ‘ imp . unconverted '  ) 

************************************************************** 

*  IF  RANDOM  MODE  AMPLITUDE  IS  DESIRED  SET  RANDAMP  EQUAL* 

*  TO  1 

*  IF  RANDOM  MODE  PHASE  SHIFT  IS  DESIRED  SET  RANDPHS  * 

*  EQUAL  TO  1 


randamp  =  0 
randphs  =  1 

★  ★★★★★★★★★Hr*********-*********************-******************* 

*  INPUT  NUMBER  OF  NODES  TO  BE  MODIFIED  AND  WRITE  NUMBER* 

*  TO  OUTPUT  FILE  IMP . UNCONVERTED 

***************************************************************** 

m  =  921 
write ( 15 , 98 ) m 
98  format (lx, i5 ) 

★  ★★★*★*•★****★★★★★**★*★★*****★*★■*★**********■*•**■*★******■**■***★ 

*  INPUT  THICKNESS  OF  SHELL  * 

*  ★★★★★It*****************************-********-***-*****-******** 

h  =  0.06 

************************************************************ 

*  SET  THE  VALUE  OF  PI 

**********************************Hr******************»****** 

pi  =  3.14159265359 

•fr*********************************************************** 

*  INPUT  THE  WEIGHTING  COEFFICIENTS  OF  THE  FIRST  TEN  * 

*  MODE  SHAPES  * 

**★****•**■**************  +  •*■********■*********■***********■******* 

al  =  0.01 
a2  =  0.01 
a3  =  0.01 
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a4  = 

0, 

.  01 

a5  = 

0. 

.  01 

a6  = 

0. 

.  01 

a  7  = 

0. 

.49 

a8  = 

0. 

,  64 

a9  = 

0  , 

.  81 

alO  = 

-  1.0 

*■»■**********************************»,■****************,**,** 

*  INPUT  PHASE  SHIFT  FOR  EACH  MODE  SHAPE 

************************************************************** 

psl  =  0.0 
ps2  =  0.0 
ps3  =  0.0 
ps4  =  0.0 
ps5  =  0.0 
ps6  =  0.0 
ps7  =  0.0 
ps8  =  0.0 
ps9  =  0.0 
pslO  =  0.0 

★  ★★★★★★*★★★★★**★*★**★■***■***■***★★***■**★*★★■**★★★★★**★★★★★*★■*■*★ 

*  INPUT  SEED  NUMBER  FOR  RANDOM  NUMBER  GENERATOR 


seed  =  1.0 

********************************************************************* 

*  CALL  SUBROUTINE  SRAND  TO  GENERATE  SEED  NUMBER  FOR  * 

*  RANDOM  NUMBER  GENERATOR  SUBROUTINE  * 

******'***************-**•***'*********************-**-*'********** 

call  srand(seed) 

*  CALL  SUBROUTINE  RAND  TO  GENERATE  RANDOM  PHASE  SHIFT  * 

*  FOR  EACH  MODE  SHAPE  * 

***★*********************************,.****.**.*.***********.**** 


:  ( randphs 

. eq. 1 ) 

then 

psl  =  2. 

*  pi  * 

rand ( ) 

ps2  =  2. 

*  pi  * 

rand { ) 

ps3  =  2. 

*  pi  * 

rand ( ) 

ps4  =  2. 

*  pi  * 

rand { ) 

ps5  =  2. 

*  pi  * 

rand ( ) 

ps6  =  2. 

*  pi  * 

rand ( ) 

ps7  =  2. 

*  pi  * 

rand ( ) 

ps8  =  2. 

*  pi  * 

rand ( ) 

ps9  =  2. 

*  pi  * 

rand ( ) 

pslO  =  2. 

*  pi  * 

'  rand ( 
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endif 


************************************************************ 


*  CALL  SUBROUTINE  RAND  TO  GENERATE  RANDOM  WEIGHTING  OF  * 

*  MODE  SHAPES  * 

ft*******************************-*******************'********* 


i f  ( randamp 


al  = 

al 

* 

a2  = 

a2 

* 

a3  = 

a3 

★ 

a4  = 

a4 

★ 

a5  = 

a5 

★ 

a6  = 

a6 

* 

a  7  = 

a7 

★ 

a8  = 

a8 

★ 

a9  = 

a9 

it 

alO  = 

=  alO 

endif 


.eq.  1)  then 
rand  ( } 
rand  ( ) 
rand ( ) 
rand ( ) 
rand  ( ) 
rand ( ) 
rand  ( ) 
rand  ( ) 
rand  ( ) 

*  rand ( ) 


*★***★★★★*★*★*★★*******■***★★★★★★*★★★★★**★★★★★★★★*•*■**★****★■*★ 

*  DO  LOOP  TO  READ  INPUT  DATA  * 

************************************************************ 


do  100  1  =  1,111 

read(14,*)  n ,  bl ,  x,  y ,  z  ,  b2 

*  CALCULATE  RADIUS  OF  NODE  POSITION 

*  (Look  carefully  at  the  axis  of  rotation  of  the  * 

*  cylinder  to  determine  the  proper  Cartesian  * 

*  coordinates  used  to  calculate  the  radius  of  the  * 

*  cylinder.  In  this  case  the  axis  of  rotation  is  the  z* 

*  axis.  * 

************************************************************ 


r= (x*x+y*y ) **0.5 


************************************************************ 

*  CALCULATE  THE  ANGLE  FROM  THE  REFERENCE  AXIS  TO  THE  * 

*  NODE  POINTS  (Look  carefully  at  the  zero  angle  * 

*  direction.  In  this  case  the  x  axis  points  in  the  * 

*  direction  of  the  charge.)  * 

************************************************************ 

if  ( (dabs (y ) . It .  0.001))  then 
if  (x.gt.0.0)  then 
a=0 . 0 
else 
a=pi 
endif 
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★  ★  ★  ★  ★ 
★ 

★ 


★  ★  ★  *  ★ 
* 

★  ★  ★  ★  * 


★  ★  ★  ★  ★ 
★ 

998 

100 

★  ★  ★  ★  ★ 
★ 

★  ★  ★  ★  ★ 


a=datan (y/x) 
if  (x. It . 0 . 0 )  then 
a=a+pi 
endif 
endif 

***********lf*l(-*1«-******»**»*****lr1r*********»*1lr****'»1r*** 

MODIFY  THE  NODE  POSITIONS  USING  THE  FIRST  TEN  MODE 
SHAPES 


* 

★  ★  * 

★ 

★  *»★*****★*★★* 

drl 

=  al 

★ 

h 

★ 

dcos(a  +  psl) 

dr2 

=  a2 

★ 

h 

•* 

dcos (2 .  *  a  +  ps2 ) 

dr3 

=  a3 

★ 

h 

★ 

dcos 1 3 .  *  a  +  ps3 ) 

dr4 

=  a4 

★ 

h 

★ 

dcos ( 4 .  *  a  +  ps4 ) 

dr5 

=  a5 

★ 

h 

★ 

dcos(5.  *  a  +  ps5) 

dr6 

=  a6 

★ 

h 

★ 

dcos(6.  *  a  +  ps6) 

dr7 

=  a  7 

★ 

h 

/ 

7 .  /  7 .  *  dcos ( 7 .  * 

a 

ps7  ) 

dr8 

=  a8 

★ 

h 

/ 

8 .  /  8 .  *  dcos ( 8 .  * 

a 

+ 

ps8 ) 

dr9 

=  a9 

★ 

h 

/ 

9.  /  9.  *  dcos (9.  * 

a 

+ 

ps9 ) 

drlO 

=  al 

0 

★ 

h 

/  10.  /  10.  *  dcos ( 

10 

* 

a  +  pslO) 

dr  = 

drl 

+ 

dr2 

+  dr3  +  dr4  +  dr5  + 

dr6 

+  dr7  +  dr8  + 

dr9  +  drlO 


★  ★*★★★*★*•*★■**★★★★★★★★★★■*★*★★★★•*★**★★★*'****★**★***•*•★**•*■ 

CALCULATE  THE  CARTESIAN  COORDINATES  OF  THE  NODES  * 

Dr*****-****-******************************************** 

x=x-dr*dcos (a ) 
y=y-dr*dsin (a ) 

WRITE  DATA  TO  TEMPORARY  DATA  FILE  IMP . UNCONVERTED 

*★***★★★■*★★★★*★★★■*★★■*★*★★■*★★★*•**★■*★★*★*■*★★•*■**■*★★*•★**★★ 

write ( 15 , 998 )  n, bl , x,y , z, b2 
format (i8,f5.0,3e20.13,f5.0) 
continue 

CLOSE  INPUT  AND  OUTPUT  FILES  * 


close  ( 14 ) 
close (15) 
stop 
end 
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program  convert imp 

★  ★*★*■*•*•**★★***-*■*■****-***■*■*•*•★*■****■•*★*****★****■*****'****★•*-***** 

*  This  program  written  by  Donald  T.  Hooker  II  on  3/15/93  * 

*  Revised  3/30/93  * 

*  This  program  converts  numerical  data 
from" imp . unconverted" * 

*  into  ascii  format  for  use  in  modifying  ingrido  for  modal* 

*  imperfections  * 

****************•*-#*********«**-*********-***-**-****«**«-•-**'**'*-** 

*  VARIABLE  DECLARATION 

*  f*  -  data  to  be  converted  to  ascii  format  * 

*  m  -  number  of  nodes  * 

*********************************************-*-*******-******* 

character*13  fl 
character*!  f2,f3,fz 
character*18  f4 
character*l  f5,f6 
character*18  f7 
character*l  f8(f9 
character*18  flO 
character*5  fll 

open ( 14 , f ile= ' imp. unconverted ' ) 
open (15, f ile= ' imp . dat ' ) 
f z= ' 0  ' 


READ  NUMBER  OF  NODES  TO  BE  MODIFIED 


read ( 14 , * )  m 


DO  LOOP  TO  READ  DATA  AND  CONVERT  IT  TO  ASCII  FORMAT  * 


do  100  i=l,m 

read (14, 998) fl, f 2, f 3 , f4, f 5, f 6, f7, f8, f 9, f 10, f 11 
998  format (al3,al,al,al8,al,al,al8,al,al,al8,a5) 

write  (15, 998)  fl,  f3,  fz,  f4,  f  6,  fz,  f7,  f  9,  fz,  flO,  fll 
100  continue 


*  CLOSE  INPUT  AND  OUTPUT  FILES  * 

•fr*********************************************************** 

close ( 14 ) 
close (15) 
stop 
end 
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APPENDIX  B 


INGRID  INPUT  FILE  FOR  TWO  DIMENSIONAL  INFINITE  CYLINDER 

MODEL 


INFINITE  CYLINDER  MODEL 

dn3d  vec  term  5.0e-3  plti  1.0e-5  prti  1000. 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.060  endmat 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  2 . 5e-8  1.0  2.5e-8 


plan  2 

00  -.003  0  0-1  0.001  symm 

00  .003  0  0  1  0.001  symm 

start 

-1  6  -11  ; 

-1  6  -11  ; 

1  2  ; 

-1.  0.  1. 

-1.  0.  1. 

-.003  .003 

all033036.0 

pri  -1-3  ;  -1-3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 
end 


89 


APPENDIX  C 


MATERIAL  MODEL  DETAILS 

For  the  steels  used  in  the  numerical  models,  a  bilinear 
stress-strain  type  model  was  used.  The  material  model  used 
was  a  "Kinematic/Isotropic  Elastic  Plastic"  model  (Hallquist 
1990).  Figure  82  shows  a  typical  stress-strain  diagram  for 
the  materials  used.  A  hardening  parameter  of  0,  giving  pure 
kinematic  hardening,  was  used  for  all  of  the  steels.  Strain 
rate  hardening  was  not  used  for  the  materials  in  this  study. 

Typical  properties  of  the  materials  used  in  the 
numerical  models  are: 


Mild  Steel  HY-100 

Density(p)  (lbm/ft3)  :  490  490 

Poission's  Ratio(v)  :  0.3  0.3 

Yield  Stress(Oy)  (psi)  :  3.2x10^  1.08xl05 

Young's  Modulus(E)  (psi)  :  2.9xl07  2.9xl07 

Hardening  Modulus (Etan!  (psi)  :  5 . lxl 04  5 . 02x1 04 


For  the  deformations  encountered  in  this  study  the  materials 
did  not  approach  their  fracture  stress. 
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Stress  (ksi) 


Figure  82.  Typical  Bilinear  Stress-Strain  Diagram 
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APPENDIX  D 


INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  STIFFENED  INFINITE 

CYLINDER  MODEL 


STIFFENED  INFINITE  CYLINDER  MODEL 

dn3d  vec  term  5.0e-3  plti  1.0e-5  prti  1000. 

c 

c  material  definitions 

c 

c  shell  material  (mild  steel) 

c 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.060  endmat 
c 

c  stiffener  material  (mild  steel) 

c 

mat  2  type  3  e  2 . 9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.12  endmat 

c 

c  load  curve  definitions 
c 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  1.0e-6  1.0  1.0e-6 

c 

c  symmetry  planes 
c 

plan  2 

00-6.  0  0-1  0.001  symm 

00  6.  0  0  1  0.001  symm 

c 

c  construct  shell 

c 

start 

-1  11  -21  ; 

-1  11  -21  ; 

1  6  11  ; 
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-1.  0.  1. 

-1.  0.  1. 

-6.  0.  6. 
a  1  1  0  3  3  0  3  6.0 

pri  -1  -3  ;  -1  -3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 
c 

c  construct  stiffener 

c 

start 

1  2  12  22  23  ; 

1  2  12  22  23  ; 

-1  ; 

-1-1011 

-1-1011 

0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

allO  550  3  6.0 

a  2  2  0  440  3  5.0 

d  2  2  0  4  4  0 

mate  2 

end 

end 
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APPENDIX  E 


INGRID  INPUT  FILE  FOR  THREE  DIMENSIONAL  STIFFENED  FINITE 

CYLINDER  MODEL 


STIFFENED  FINITE  CYLINDER 

dn3d  vec  term  5.0e-3  plti  1.0e-5  prti  1000. 

c  material  definitions 

c  shell  material  (mild  steel) 

mat  1  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.060  endmat 
c 

c  stiffener  material  (mild  steel) 

c 

mat  2  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.1e+4  sigy  3.2e+4  shell  quad  5  thick  0.12  endmat 

c  endplate  material  (HY-100  steel) 

mat  3  type  3  e  2.9e+7  pr  0.3  ro  7.356e-4 

etan  5.02e+4  sigy  1.08e+5  shell  quad  5  thick  0.25 
endmat 

c 

c  load  curve  definitions 
c 

led  1  2  0.0  0.0  1.0  0.0 
led  2  2  0.0  1.0e-6  1.0  1.0e-6 

c 

c  symmetry  planes 

c 

plan  1 

00  0.0  0  0-1  0.001  symm 

c 

c  construct  shell 

c 

start 

-1  6  -11  ; 
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-1  6  -11  ; 

1  16  ; 

-1.  0.  1. 

-1.  0.  1. 

0.0  18.0 
a  1  1  0  3  3  0  3  6,0 

pri  -1  -3  ;  -1  -3  ;  ;  1  -1.0  0.  0.  0. 
mate  1 
end 

c 

c  construct  stiffener 

c 

start 

1271213; 

1  2  7  12  13  ; 

-1  ; 

-1-1011 

-1-1011 

6.0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

allO  550  3  6.0 

a  220  440  3  5.0 

d  2  2  0  4  4  0 

mate  2 

end 

c 

c  construct  endplate 

c  surface  definitions 

c 

sd  1  cyli  000  001  6.0 

sd  2  cyli  000  001  [6. 0*3/5] 

start 

1  5  10  15  19  ; 

1  5  10  15  19  ; 

-1  ; 

-1-1011 

-1-1011 

18.0 

prill  551  1  -1.0  0.0  0.0  0.0 

di  1  2  0  4  5  ;  1  2  0  4  5  ;  ; 

sfvi  -2  -4  ;  -2  -4  ;  ;  sd  2 

sfi  -1  -5  ;  -1  -5  ;  ;  sdl 

mate  3 

end 

end 
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